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SUMMARY.   In  continuation  of  the  authors'  paper  (I966), 
here  a  class  of  rank  order  tests  for  the  homogeneity  of 
several  dispersion  matrices  is  proposed  and  their 
properties  are  studied. 

1 .   INTRODUCTION 

The  rank  permutation  technique  in  multivariate 
analysis  is  of  quite  recent  origin,  and  is  mostly  confined 
to  the  problems  of  location  and  scale,  and  of  association 
patterns  (see  Puri  and  Sen  (I965)  and  Sen  (I966)  for 
references).   But  the  authors  are  not  aware  of  any 
nonparametric  test  for  homogeneity  of  dispersion 


matrices.   The  object  of  the  present  Investigation  is  to 
develop  some  nonparametric  procedures  for  the  problem 
of  testing  homogeneity  of  several  dispersion  matrices 
through  the  use  of  a  class  of  rank  order  statistics  which 
are  extensions  of  Chernoff  and  Savage  (1958)  type  of 
statistics . 

In  the  parametric  theory,  the  problem  of  testing 
homogeneity  of  dispersion  matrices  has  nothing  to  do  with 
the  homogeneity  of  location  parameters.   However,  in  the 
nonparametric  case  the  formulation  of  the  hypothesis  of 
invariance  deems  the  identity  of  distributions,  and  this, 
in  turn  requires  the  homogeneity  of  location  parameters. 
Thus  we  consider  the  following  problems:   (i)  Tests  for 
the  homogeneity  of  dispersion  matrices  assuming  the  identity 
of  location  parameters.   Some  permutationally  distribution 
free  tests  for  this  problem  are  considered  here  and  their 
large  sample  properties  are  studied;  (ii)  Tests  for  the 
homogeneity  of  dispersion  matrices  without  assuming  the 
identity  of  location  parameters.   Here  the  basic  permutation 
argument  is  not  tenable,  and  only  large  sample  properties 
are  considered.   Finally  (iii)  the  problem  of  testing 
simultaneously  the  homogeneity  of  location  parameters  as 
well  as  the  dispersion  matrices  is  briefly  discussed. 

In  the  parametric  case,  tests  for  the  homogeneity 
of  dispersion  matrices  are  based  on  the  sample  covariance 


matrices  [see  Anderson  (I958,  Chapter  10)].   However,  in 
the  nonparametrlc  set  up,  the  sample  covariance  matrix 
Is  not  usually  looked  upon  as  a  suitable  statistic  (vN^rv  Xvc  a  ) 
because  of  its  sensitiveness  to  outlying  observations,  of 
the  necessity  of  the  existence  of  the  moments  of  the 
parent  distribution,  and  of  the  fact  that  it  may  not 
retain  its  sufficiency  or  optlmality  for  non-normal  distri- 
butions.  Due  to  these  reasons,  we  shall  formulate  a  class 
of  regular  functionals  which  seems  desirable  for  nonpara- 
metrlc procedures.   The  corresponding  rank  order  estimates 
are  then  used  to  formulate  the  proposed  tests. 


2.   PRELIMINARY  NOTIONS 


Let  X^^^=  {x[^^  ,  .  .  .  ,X^^h  ,    a  =  l,...,n,  be  n, 
^a  ^  la  '    '  pa  ' '       '    '  k     k 

independent  and  identically  distributed  (vector  valued) 

random  variables  (i-i-d-r-v)  having  a  p-variate  continuous 

cumulative  distribution  function  (cdf)  F^'^^(x),  for 

k  =  l,...,c(>  2);  all  these  c  samples  being  assumed  to  be 

c 
mutually  independent.   Denote  N  =  >   n,  and  assume  that 

k-1   ^ 
for  all  N,  the  Inequalities 

0  <  X^  <  aJj^^  =  Hj^/N  <  1  -  A^  <  1  ,   for  k  =  l,...,c   (2.1) 


hold  for  some  fixed  A^  <  l/c . 

o  = 


Let  F)'^l(x)  and  P).  .Jx,y)  be  the  marginal  cdfs  of 

x!^^  and  (x!^^,X^.^h  respectively  for  1  <  j  =  l,...,p,  and 

k  =  l,...,c.   Let-'?'  be  the  class  of  all  continuous 

-  p 

p-varlate  cdfs  for  p  >  1.   Then,  we  assume  that 


,(k) 


,(k) 


(k) 


^V  "^  -^p  =^  F(i)(^)'^  ^1  and  Fjj;j)(x,y)  <^  .l^'^  •  (2-2) 


Denote 


H(i)(^)  =|Z4^^  ^(1)(^) 


(2. 3) 


Thus  H/.^(x)  is  the  combined  population  cdf  of  the  ith 
variate  for  i  =  l,,..,p.   Let  J/.n(u)  for  i  =  l,...,p  be 
some  absolutely  continuous  function  of  u,  defined  in  the 
open  interval  (0,1)  and  suppose  that  J,.,{u)  is  normalized 
in  the  following  manner: 


0 


J/ . ^ (u)  du  =  0 

(l) 


0 


J,^x(u)  du 


=  1;  i  =  l,...,p   (2.4) 


Later  on  we  shall  impose  certain  regularity  conditions 

on  J/.v(u).   Let  us  define 
(i) 


.(k) 
ia 


(1' 


H,.,(x!^^ 
( 1 ) ^  la  ' 


i  =  1, . . . ,p  ;  k  =  1, . . . ,c.  (2.5) 


and  we  shall  call  Y|^'  as  the  generalized  grade  of 

(k) 
■^la    ^^  shall  employ  these  generalized  grades  to 

formulate  the  desired  class  of  tests. 

Let  us  now  define 


0^j(H)  -  GOV 


(J(i)[H(,)(x)],  J^.)[H(.^(y)]) 


H{^) 


OO  oo 


J(i)[H(.)(x)J  J(.)  [H(.)(y)J  dH(.^.^(:<,y)   (2.6) 


-00  -00 


where 


H/  .  . ^ (x , y ) 
(1,  j)  '^ 


=  V 


k=l 


(k)   (k)   ,   ^> 


for  i  /  J  =  1,. . .,p   (2.7) 


and  H,^^(x)  is  as  defined  in  (2.3).   Also  denote 


e[f  -  cov 


(J(i)[H(,)U)],  J(j)[«(J)(^)]^ 


F^^^(^) 


00  00 


J 


(i) 


H(,)(.) 


-00  -00 


I  ■\       H,  .,  (y) 
(  J  )  L  (  J  )    _ 


(2.8) 


where 


M. 


(k) 
N,i 


00 


(i) 


H(,){x) 


dFJJj(x); 


i  =  1,. . . ,p  ;  k=l,...,c  (2.9) 


-00 


Finally  let 


(Sf  =  {{6    AE)))^  and  ©^^^  =  (  ( ej  ^  )  )  i  1-1     n^ 


k=l,...,c.  (2.10) 


Thus,  @^    is  the  dispersion  matrix  of  the  generalized 
grades  of  the  kth  sample,  k  =  l,...,c.   The  purpose  of 
introducing  the  dispersion  matrices  in  (2.10)  is  to  work 
with  a  class  of  dispersion  matrices  which  are  invariant 
under  certain  transformations  on  the  variables,  are  less 
sensitive  to  outlying  observations,  and  at  the  same  time, 
are  reasonably  informative  measures.   In  particular,  if 
J,.s(u)  is  monotone  in  u  ( 0  <  u  <  1 )  for  all  i  =  1, . . . ,p, 
(H/   ,  k  =  l,...,c  defined  in  (2.10)  will  be  invariant 
under  monotone  transformations  of  the  co-ordinate  variables. 

Now  since  the  cdfs  F^  ;...,F^'^^  are  all  unknown, 
so  are  also  the  (c+1)  matrices  defined  in  (2.10).   To 
estimate  these  we  follow  Von  Mises  (19^7)  approach,  and 
define  for  any  q  <  p 


<  (x-j_,  .  .  .,x^)  ,  a  =  l,...,n^]  (2.11) 


%il,...,iq)^^l'---'\)  ^I^'n^^  ^N(V...,i  )[-!'- ■■'\U2.12) 


for  all  in  ?^  •••/^  ig  ^  1>---jP;  Q.  =   1.  •••>?•   Again,  let 

AV^,...,1    )(^i---^q)  -  PM'^a  ^  ^J  '    J  =  l>---^5  (2.13) 


"(1-^,...,!  J^'^l'-'-'^q^ 


(k)  ^(k) 


^^N    ^U;,...,!  )(^l"-"^q)  (2.14) 


'or  all  i,  ^  ...  ;^  i„  =  !»•••>?;  Q  =  1?  •••>?■ 


Now  in  actual  practice,  we  have  a  sequence  of  functions 
J^/.v(u)  which  converges  to  J,.^(u)  as  N  — >  oo  for  all 
0  <  u  <  1,  and  i  -   l,...,p.   Thus  we  may  estimate  0.  .(H) 
defined  in  (2.6)  by 

oo    00 

'  f  J 


T     - 
N,ij 


N(i) 


"N(i)(^^ 


'N(j) 


[H^(.)(y)]  dHj,(,^.)(x,y) 


-oo  -00 


C  k 


^N  ^S^  JN(i)[%(i)^^ia  )J  h{j) 


%(jV^Ja^)J  (^-^5^ 


for  all  i,J  =  l,...,p.   Since  T^  ^  ,•  =  T»,  . .  ,  the  number 
of  statistics  T^  .  .  is  actually  p(p+l)/2.   Now  let  us 


denote 


^N  =  K^^.  ,  i  <  j  =  1,...,P} 


(2.16) 


00 

^N,i  =  J    '^N(i)L"N(i) 
-oo 


J»'' j^)  Hj^j/  ^  ^  (x 


dF^'j^j)(x)  ;   i  =  l,...,p;  (2.1?) 


■N,lj 


oo  oo 


-OO  -CXD 


'N(i) 


8 


HN(i)(^) 


^N(J)  L"N(j)(y\ 


dF(^l,j)(.,y) 


n 


k 


k  a=l    ^ 


_%(l)(^ia  y  "^N(j)  L^N(j)(^ja  ^J 


m(k)   (k) 
"   N,l  ■'n,j 


(2.14) 


and,  finally 


T 


(k; 

N 


(k) 


IN,! 


1  < 


J  =  1,. . .,p];   k  =  1, . 


,c 


(2.19) 


Then,  our  proposed  tests  are  based  on  some  functions 
(see  (3-24) )  of  tJ 


(k)  * 

'^    '    ,    k  =  1,  .  .  .  ,c  and  T, 


It  may  be  noted  that  the  class  of  dispersion  matrices 
considered  above  contains  many  well  known  nonparametric 
measures  as  special  cases.   For  example,  if  we  let 


Jj^(^^  (a/(N+l)) 


12/(N^-l)^/2 


(a  - 


N+1 


)  ,  a  =  1, . . . ,N; 


i  =  l,...,p;then  T^  reduces  to  the  Spearman's  rank  correlations 
[cf.  Kendall  (1955)],  while  for  J,^Ju)  =  v'12"  (u  -  l/2), 
(gT  reduces  to  the  grade  correlations  matrix.   The  nonpara- 
metric  measures  by  Blomquist  ( I95O) ^  Bhuchongkul  (1964) 
among  some  others  belong  to  this  class. 

We  now  consider  the  statistical  formulations  of 
the  problem  of  testing  the  homogeneity  of  dispersion  matrices. 


First  we  consider  equality  of  covariance  matrices,  assuming 
the  identity  of  locations.   In  such  a  case  our  null 
hypothesis  is 

jj(l).   f(1)(x)  =  ...  ^  F^^\x)  (2.20) 

(which  implies  that(gf''"^  =  ...  =(g)^'^^),  against  the  set  of 
alternatives  that  @^-^' ,  ...  J^^^    are  not  all  identical. 
We  may  note  that  H^  '    will  include  test  for  identity  of 
scales  and  of  association  patterns  of  the  c  cdfs.   Next 
we  consider  the  equality  of  dispersion  matrices  without 
assuming  the  identity  of  locations.   In  such  a  case, 
denoting 

F^^^(x)  =  f[^Ux    -   5^^^)  ,  k  =  l,...,c      (2.21) 

Where  6^^^  =  (5^^^ . . . ,5^^^ )  ;  k  -  1, . . . ,c  are  c  real 
p-vectors,  we  test 


H^^^  =  F^^^x)  =  ...  =  Fi^^(x)  (2.22) 

o      o   ~  o   ~ 


against  similar  type  of  alternatives  involving 

(B^''"  ,.../H^'^^,  where  of  course,  in  the  definition  of 

(j^<k)  ,g  .  f(1)^.  _^p(c)are  replaced  by  fI^\  . .  .  ,fI^'' 

respectively.   Finally,  we  consider  the  hypothesis  that 
both  locations  and  dispersion  matrices  are  the  same;  that 
is,  we  test 
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^(3).   f(1)(^)  =  ...  =  F^^^x)  (2.23) 


where  F^   (x),  k  =  l,...,c  is  given  by  (2.21),  against 
the  alternatives  that  at  least  one  of  the  following  two 
inequalities  is  not  true. 

6^1)  =  ...  =  5(^)  ;  (S^l^  ^  ...  =(H)(^)  .      (2.24) 


ru  ■  O/ 


3.   PERMUTATION  TESTS  FOR  H^"^^ 
-^  o 

Let  us  pool  the  N  =  y   n,  observations  X^   , 

k=l 

a  =  l,...,n,  ;  k  ~   l,...,c  into  a  combined  sample,  and 
denote  the  sample  point  by  p  x  N  matrix 

^N  =  (^r^'""Sn'^')  (5.1) 

c 

and  the  sample  space  by  'K..      Under  H^   ,  the  Joint  distri- 
bution of  Z„  remains  invariant  under  the  finite  group  of 
NJ  permutations  of  the  N  columns  of  Z  .   Hence,  conditionally 
given  Z^,  all  the  N!  sample  points  generated  by  these  NI 

T 

permutations  of  the  columns  of  Z,,  are  equiprobable,  each 
having  the  conditional  probability  l/N! .   Let  us  denote 
this  conditional  probability  measure  (defined  over  the  NI 
equiprobable  realizations)  by  X^.   Now,  ranking  the  N 
elements  in  each  row  of  Z,,  in  increasing  order  of 


11 


magnitude,  we  obtain  a  p  x  N  matrix 


^N  = 


R 


(1) 
11 


\r(^) 


pl 


R 


(1) 
In. 


R 


(1) 
■pn. 


R 


(c) 


(3-2) 


where  by  virtue  of  continuity  of  the  cdfs,  the  possibility 
of  ties  may  be  ignored  in  probability.   For  every 
i(=  l,...,p),  replacing  the  ranks  a(=  1,...,N)  in  the  ith 
line  of  R-j  by  a  set  of  general  scores 


-■"N  a  ^  "^Nf  i)  ^^'^^^ '  ^   ^   l>---^NJj  we  get  the  correspondin 


g 


score  matrix 


En  = 


/^(l) 


N,R 


(1) 
In. 


,(P) 

N,R 


(1) 

pn-i 


N,R  ='  \ 

'  In   \ 
c  , 


Ap) 

N,R 


(c) 


(5.3) 


(K) 
With  this  terminology,  we  can  write  T,-,    .  .  and  T.,  .  . 

^^ '  N,ij      N,ij 


defined  by  (2.l8)  and  (2.15)  as 


(3.4) 


n 


k 


(i) 


,(k) 


.(j) 


n, 


n 


,p(l<:)  ^_1_ 

^,^^      "k  d=r  N,R.''''   NjRl'" 
'  ia      Ja 


(k) 


^  V— E 
n 


(i) 


(k) 


k  a=l  N,R 

'  la 


k 


(j) 


k  a=l  N,R 


(k) 


1  <  j  =  l,...,p;  k  =  1,...,C 
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n 


■N,1J    N 


T 


k 


.(1)    Jj) 


k=l  a=l  N,R 


(k) 
ia 


'  act 


i  <  J  =  l,...,p  .  {3.5) 


Operationally,  It  will  be  much  simpler  to  work  with 
the  slightly  adjusted  statistics: 


.(k) 


n 


k 


(i: 


^(j) 


'^J   V^^S^N,RP^)   N,R^^ 
'  la      30. 


n 

^kferN,R!^L 


^-^--  -(k) 
J  a 


kx=l  N,R 


}  (3.6) 


i  i  J  =  lj--',P;  k  =  l,...,c  , 


n 


N,lj 


N-li 


k=T  a=r  N,r!^^   N,R^.^^ 


Z 


N  E^ji)  E^JH      (3.7) 


where 


E 


(1: 

N 


f5 


k=l  a=l  N,R 


(k) 
ia 


0 


'N(i; 


HN(i)(^M^\(i)(^^i=l^---.P  (5.8) 


Note  that 


and 


n/(n^-l) 


T 


(k) 
N,ij 


N,lj 


N/N-1) 


\,ij  -  ^4"^  4'^ 


(3.9) 


Now  under  the  permutational  probability  measure  /,,, 
there  are  NI  equally  likely  (conditionally)  permutations 
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of  the  columns  of  Rj,  (and  hence  of  E„).   Hence,  by  an 
adaptation  of  the  same  argument  as  in  Purl  and  Sen  (I966), 
we  have 


E(si^!  . 
N,ij 


'^N^  "  ^N,ij   ^°^  all  i  <  j  =.  l,...,p 


,(k) 
'N,lj 


Furthermore,  S^,  :.   in  (3-6)  can  also  be  expressed  as  a 


U-statlstics  in  the  scores  e'"^^,  >  ;  i  =  l,,.,,p; 
a  =  l,...,n,  as  follows: 


n, 


;(k)   _  1  ("k,-l  ^  ,  (1)       (1)    ,  (J)       (J) 

J-U  ±p         ju  jp 

(3.10) 

for  i,  j  =l,,..,p;  k=l,...,c  . 


Let  us  now  define 

-I  c         c         c  k        q        r         ,  .  ■.  ,  .  > 

N(N-l)(N-2)    ^IZ>Z:iZ5ZIZ   (E^'\k)    -    ^^       (q))     ' 

fgCj)  e(J)  ^     .     /pd')  e(I')        wp(j')  p(j')        ^ 

N  r(^^         n  r(^)  n  r(^)   "    N  r(^)       N  r(^)   "    n  r(^) 
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1   ^ 

N(N-l)    -—J 

c 

^k      ^k 

■^irr 

(E 


(1) 


,(k) 


N,R^)    N,R^) 


-  E 


(j) 

N,R 


(q: 


,  (IM        (!•)    Wp(j')       p(j')    ^     c.*      q*  (^    lO) 


for  1,  j,  i',  j'  =  1, .  .  .,p 


Then,  by  an  adaptation  of  the  results  of  Nandi  and 
Sen  (1963)  and  Sen  (I966),  It  can  be  shown  after  some 
algebraic  manipulations  that 


n,  -2  n  -2 

(  \  )    2     p   aJ".   ) 


n,    ^-^-  '  d'  I  n,  -2 


q  ,  (5.15) 


i'^fr 


1  ^  (2)(N-2.   (d)     ,   , 


if  k  ^  q,      (3.14) 


for  i,  1',  J,  j'  =  l,...,p;  k,q  =  l,...,c  . 

Thus,  for  small  samples  one  may  work  with  the  random 
variables 


l^Nfij  -  ^N,lj  '      1  ^  J  =  1,--.,P.  k  =  l,...,c}  ,  (3.15; 


IS 


and  considering  the  generalized  inverse  of  their  (permutation) 

covariance  matrix,  arrive  at  the  desired  quadratic  form 

as   suitable  test  statistics.   Under  ?^   this  statistic 

c 
will  have  only  Nl/rPn,  I  equally  likely  (not  necessarily 

all  distinct)  values,  and  hence,  an  exact  size  e(0  <  e  <  1) 
test  for  }i\^        can  be  constructed.   However,  the  expressions 
in  (3-11)  through  (3.14)  are  quite  complicated.  If  N  is 
not  too  small,  it  will  be  convenient  to  consider  the 
following  permutation  test. 

By  virtue  of  (2.4),  E^^'    defined  in  (3-8)  converges 
to  zero  as  N  — >  oo .   In  fact,  (we  shall  see  later  on) 
E^^^  -  o(N"^/^)  for  all  i  =  l,...,p.   Hence  from  (3.5) 
and  (3.7) 


*       N-1  ^ 
^N,ij  ~   N   ''N,ij 


=  o(N-l/2) 


(3.16) 


Further,  from  the  results  of  Puri  and  Sen  (196^),  under/' 


N 


n. 


-ZZE(i\^,  -li^M  =  0  (N-1/2) 
'  la 


(3.17) 


for  all  i  =  l,...,p;  k  =  l,...,c.   Consequently  under ^Jl 


n 


,(k) 
'N,ij 


k 


n 


2_ 


.(i) 


k  a=l   N,R 


(k) 
la 


r(j) 
N,R 


(k) 


=  Op(N-^) 


(3.18) 


for  all  i,  J  =:  l,...,p;  k  =  l,...,c  . 
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Now  let  us  define 


c    k 


V. 


n- 1  ^  i}.w  >         T\T  <: ^ 


(1)     ^(j) 


E^^'!,  .  e(J') 


UM'j'^^N'  -N^^-^^j^U)  -^^j,(k)  -^^  (k)  ^M 


la 


ja 


1  'a 


J  'a 


(5.19) 

for  all  1,  i',  J,  j'  =  l,...,p  . 

Setting 

r  =  (i  -  l)(2p  -  i)/2  +  j   for  i  <  j  =  1,...,P  (3.20) 


we  rewrite 


.^N,ij  '    ^    =    ^ 


l,...,p}   as  {S^^^  ,  r  =  l,...,p(p+l)/2  \^   S^^^  (3. 21) 


for  k  ==  1, . . .  ,c  , 


1^^,13    '    '    =   ' 


=  l,...,v]     as  [S*  ^  ,  r  -  l,...,p(p+l)/2}  =  S*   (3.22) 


and 


\(\^  =    ((''rs^^N^^^r,s  =  1,  .  .  .  ,p(p+l  )/2       (3-23) 


where  r  is  defined  as  in  (3-20)  with  (i,j)  replaced  by  (i',j') 

Then,  proceeding  as  in  Purl  and  Sen  (I966),  we  consider 
the  statistic 


^N 


X 


n 


k=l 


k 


^N    "  ^N 


V'(^n) 


o  (k)    * 
^N    "  ^N 


(3.24) 
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for  the  problem  of  testing  hJ^-"-^.   (Here  VJ^-'-(R  )  is  in 
the  inverse  of  Vj^(R,j),  and  is  assumed  to  be  positive 
definite).   Thus,  we  have  the  following  rule: 


(1)  (5.25) 

<  £^^g  ,  accept  hJ  "^ ' 


o 


where   £^        is   so  chosen  that    P{£>j  ^  £•>! 


H<^'i  ^  e. 


where  0  <  e  <  1  is  the  preassigned  level  of  significance 
of  our  test.   In  small  smaples,  we  can  evaluate  £^ 
by  referring  to  the  exact  permutation  distribution  of  Z^ 
generated  by  the  NI  equally  likely  permutations  of  the 
columns  of  {3-3),   while  for  large  samples,  we  shall  show 
in  the  next  section  that 

%,£  ~^  '^  e,(c-l)p(p+l)/2  (3.26) 


2 
e 

with  r  degrees  of  freedom. 


where  'X  ^  r.  ^^  ^^^   lOO(l-e)  y,  point  of  a  chldistribution 
e ,  r 


k.      ASYMPTOTIC  PERMUATION  DISTRIBUTION  OF  Z^^ 
As  in  Section  3>  we  define 
eJ[,^^  =  J^(^)(a/N);   a  ^  1,...,M;   i-l,...,p,   (^.1) 
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and  as  In  Chernoff  Savage  (I958),  we  extend  the  domain 
of  J-j/.<(u)  to  (0,1)  by  letting  Jja/ •  n  be  constant  on 
Ta/N,  (a+l)/N).   Furthermore,  we  shall  make  the  following 
assumptions* : 


Assumption  I.  J,.k{\i)    -   lim   J^/.^fu)  exists  for 

(^)      N— >oo   ^(^^ 


0  <  u  <  1  and  is  not  a  constant  (for  i  =  l,...,p) 


Assumption  II. 


{'mi: 


%(i)(^ 


J 


N(J) 


^N(i)''%(j) 


-  J 


(i) 


%(i)(^ 


(j) 


HN(j)(y: 


(k) 


Vj)^^^J1^^(i;j)(^'^^ 


=  0p(N-l/2) 


where 


i,j  =  l,...,p;  k  =  l,...,c    (4.2) 


-N(i) 


0  <  Hjj^^(x)  <  1  j  ;  i  =  1,  .  ..,p 


Assumption  III  J^r  •  )    =  o(N"^/  ) 


(4.3) 


As  sumption  IV. 


jj[j(u) 


Su 


J /.  Au] 


V    "(i) 


(4.4 


<  K[u(l  -  u) ] 


ar 


*The  assumptions  1  to  4  are  similar  to  those  of  Bhuchangkul 
(1964).   These  assumptions  are  somewhat  more  restrictive 
as  compared  to  the  ones  in  Purl  and  Sen  (I966).   However, 
this  restriction  may  be  ascribed  to  the  fact  that  here  we 
are  dealing  with  rank  order  covariance  statistics,  while  in 
the  earlier  paper,  we  dealt  only  with  linear  rank  order 
statistics. 
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for  r  =  0,  1,  2,  where  0  <  a^  <  l/8;  a,  =  -  1  and  a^  =  -  2, 
Now  let  us  denote 


^Nij 


00      00 


(1) 


H(,,(. 


J 


(j) 


"(jV^^.. 


-co  -00 


dF^^J  .^(x,y);  k=l,...,c; 


00      00 

-oo  -oo 


H(,)(. 


(j 


"(J)^^M    ''"(l,j)(^'y^       (^-^^ 


ij,i' j' 


00      00      oo      00 


-oo  -oo  -oo  -oo 


(i) 


H(,)(x.) 


(j) 


"(J)(^j)j 


"^(IM 


"(l')(^l' 


•^(j') 


H(..)(x.,) 


dH 


(l,j,i',j')^''l'''j'''i"''j'^ 


*         * 


"^H^^i'l'      fo^   Ij  J^i'>J'=l.  •  •  •»?•    (4.7) 


Adopting  the  suffixing  system  of  (3-20),  we  rewrite 

V. .  . , . ,  as  V   for  r, s  -   1, . . . ,p(p+l)/2  and  denote 
1 J  >  1  J        I'S 


V  -   ((Vrs^^,s  =  l,...,p(p+l)/2  • 


(4.8) 


Assumption  V.   v  defined  by  (4.8)  is  positive  definite. 
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Then,  we  have  the  following  Theorems. 

THEOREM  4.1.   Under  assumptions  I  tjo  V,  Vj^(R^),  defined  In 
(3.23),  Is  equivalent  In  probability  to  v  ,  defined  in 
(4.8),  as_  N — >oo  uniformly  In  A..   ,...,A^   and  all 

T^(l)  T7(C  )^  ^-T    . 

F^   ,...,F^  "^   -^p 

The  proof  follows  precisely  along  the  same  line  as 
In  Theorem  4.2  of  Purl  and  Sen  (I966)  with  straight  forward 
extensions,  and  hence  is  omitted. 

THEOREM  4.2.   Under  ass\;imptlons  I  to  V,  the  permutation 
distribution  of  £,.  (defined  by  ( 3  •  24 )  )  asymptotically 
reduces  to  a  chlsquare  distribution  with  ( c-1 )p(p+l )/2 
degrees  of  freedom. 

Proof.   By  virtue  of  (3.18),  in  £^  we  can  replace 

n 

oCk)   ,    „(k)   _  _1_  t;^ —   (i)      ( j)  .  <  .  _  T       „. 

^N,ij  °y  ^N,1J  -  n,  -^^  ^N,R.   ^N,R.   ^  °^  all  1  <  j  -  1,  .  .  .  ,p, 

k  =  l,...,c.   We  shall  first  prove  that  under/'' 

iN^/^(Z^^|j  -  T*^.j),  1  <  j  =  l,...,p;    k  =  l,...,cj  has 
asymptotically  a  (c-1 )p(p+l )/2  variate  normal  distribution. 

For  this,  it  suffices  to  show  that  for  any  non-null  real 

c 
5  =  (5,,...,  6  ),  such  that  >   n,  5,  =  0, 

~        1       '   C   '  ^—yr-        k  k       ' 

In  /   >   6^  Z^  j   1  <  j  =  l,...,p]  has  asymptotically 
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(under  '^)  a  multivariate  normal  distribution  of  rank 
p(p+l)/2.   Now  let  US  denote 


.(k) 


N    5,  Z'  .  . 
k^  k   N,ij 


^^ — 

k^ 


^   -^k  &■  N,RP^)   N,R^^ 
'  la      ja 


c 
k^ 


k   ^N,r  -  %,r 


(4.9) 


where  r 
and  define 


(l-l)(2p-l)/2 


+  j   for  i  <  j  =  1,. . .,p  (cf.  3.20: 


N,r,a 


n,r!^)  n,r(.^) 


(l-l)(2p-i)/2 


+  J  ,  (4.10) 


then,  we  have  N  vectors 


a 


(k)'  _  ,Jk) 


Na 


U 


1^1,1,0' 


.,(k) 
■N,p(p+l)/2,a 


,cc 


)    for  a  =  1,  .  .  .  ,n,  ; 


k  =  1, 


(4.11) 


and  under  /I,  all  possible  NI  permutations  of  these  N 
vectors  among  themselves  are  (conditionally)  equally  likely, 
each  permutation  having  common  permutational  probability  l/N! 
Thus  W^  =  [%,r'-      ^   "  l,...,p(p+l)/2  j  is  a  p(p+l)/2 
vector  where  each  element  is  a  linear  permutation  statistic 
of  the  type  considered  by  Wald-Wolfowitz-Noether-Hoeffding- 


Hajek  among  others.   Thus  by  a  direct  application  of 
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Theorem  7.I  of  Hajek  (I961),  it  follows  that  W^^  has 
asymptotically  a  p(p+l)/2  varlate  normal  distribution. 
The  rest  of  the  proof  follows  by  straight  forward  routine 
analysis.   Hence  the  theorem. 

By  virtue  of  Theorem  4.2,  Z^   ^   defined  in  (3.25) 
reduces  to 'X-  t      i]    (    +-\  ) /o   which  asserts  the  truth  of 
(5.26). 


5.   ASYMPTOTIC  NORMALITY  OF  S^^^ ,  k  =  l,...,c  FOR 
ARBITRARY  F^ ^ ^ , . . . ,F^^ ^ 


The  main  theorem  of  this  section  is  the  following. 
Theorem  3.1.   Under  the  assumptions  I  to  V  of  Section  4, 


the  random  variables  N 


1/2 


U(^)      u(^)   -U^^V^^'^   T^^)-  n(^). 
^N,ij  "  ^N,ij   %,i^^N,j'  ^N,i  ^^N,i' 


i  ^  J  =  1, 


where  T„  •  and  n.\  .  .  are  given  by  (2.16), 
_  N,i  N,ij  s £ — ^     '' 

(2.8)  and  (4.5)  respectively,  and  U^^|j  =  T^j^jj  +  T^^j  T^^^J. 

have  asymptotically  a  cp(p+3)/2-variate  normal  distribution 
with  null  mean  vector  and  finite  covariance  matri:x. 


The  proof  of  this  theorem  is  given  in  the  Appendix  A. 

We  may  remark  that  the  asymptotic  normality  of  the 
random  variables  N''"/^ ( tJj'^ j  -  n^^|  ;  i  =  l,...,p;  k  =  l,...,c) 
follows  under  slightly  weaker  conditions  [cf.  Puri  and  Sen 
(1966,  Theorem  5.I)].   On  the  other  hand  the  assumptions 
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I  to  V  are  needed  for  the  asymptotic  normality  of 


N 


1/2 


U 


(k) 


(k) 


(k)  ,,(k)- 


N,ij   ^^N,lj   ^^N,l  ^N,j. 


>       -'-   _   J   ~   -'-}•••>?>  rC— -I-  ;  •  •  •  ^  C  • 


We  shall  now  obtain  the  limiting  distributions  of 
£^  given  by  (3-24)  under  the  null  hypothesis  H^^  and 
under  a  sequence  of  specific  dispersion  alternatives  to 
be  stated  below.   In  what  follows,  we  assume  that  for  all 


A 


(k) 

N 


(k) 


exists  and  is  positive. 


k  -  1, . . . , c ,  11m 

N — >oo 

Now  under  H^   ,  it  follows  from  the  results  of  Puri 

and  Sen  (I966)  that  I  T^^|  I  -  0  (N"''-/^)  for  all  i  =  l,...,p 

and  k  =  l,...,c;  and  following  routine  but  lengthy 

computations  one  can  show  that  under  H^   ,  asymptotically 


N  cov 


(k)    ^(q) 
■^N^rs 


i^N,i 


l^)    -  ^i  ^j,rs 


(5.1) 


for  all  i,  j,  r,  s  =  l,...,p;  k,  q  =  l,...,c,  where  6 


kq 


is  the  usual  Kronecker  delta  and  v. .    is  defined  as  in 

ijrs 

(4.7)  with  the  further  simplification  that  f{^}    =  ^(i)' 
^(ilj,r,s)  "  "(i,J,r,s)  ^°^  ^11  ^  -  l,---,c;    i,  j,  r,  s  = 
l,...,p.   Consequently,  using  Theorems  4.1,  5.I  and  some 
routine  analysis,  it  can  be  shown  that  under  H^   ,  £,, 
has  asymptotically  (unconditionally)  a  chi-square  distribution 
with  (c-l)p(p+l)/2  degrees  of  freedom. 

We  now  formulate  a  class  of  dispersion  alternatives. 
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For  the  scalar  alternatives  we  let  (as  in  Purl  and  Sen 
(1966)) 

FJ^j(x)  =  F(i){[x  -  Hilsj^^;  1  =  i,...,p;  k=l,...,c  (5.2) 


where  under  H^   ,  all  &:    are  equal  to  unity.   Next,  as 

in  Sibuya  (1959),  we  consider  bivariate  dependence  functions 

O:.  defined  as 


O 


(k)  _  ^(k)   ,   ,  .  /  T,(k),  .  ^(k) 


ij 


Fjj;.^(x,y)  /  F|j;(^)  Fij)(y)        (5. 5) 


for  all  i  ;^  j  =  l,..,,p;  k  =  l,...,c.   Clearly,  O-  • 

is  a  function  of  F/ • I  and  F)  .1 ,    and  will  be  equal  to 

unity  if  X:  ^and  X\^  are  statistically  independent. 
'^     la      ja 

Various  properties  of  such  dependence  functions  have  been 

studied  by  Sibuya  (1959)  and  others,  and  in  nonparametric 

field,  they  appear  to  have  some  advantages  over  the  other 

measures.   We  denote  by  O.  .,  the  same  function  (i.e.  {5-3)) 

for  the  bivariate  cdf  F, .  .\(x,y),  and  relate 

\  1 J  J  / 


^W   "  ^i  i  +  ^iV   ^°^  i  ^  j  =  l,...,p;  k-l,...,c  (5.4) 


(5.2),  (5.3)  and  (5.4)  together  constitute  our  desired 
class  of  alternatives.   Now  we  can  rewrite  (2.8)  as 
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OO  00 


0 


(k)  _ 


(k) 


-J  J    L^l"j>'=<- 


-OO  -00 


FJ^jU)  FS5J(y) 


(1) 


H(.)(x) 


J 


(j) 


,"(J)(^' 


clH 


(.^(^)  dH(j)(y)  (5.5) 


for  all  1,  j  =  l,...,p;  k  -  l,...,c.  From  {5-3),  (5-^), 
we  note  that  0^  .  measures  In  a  sense  some  functional  of 
the  dependence  function  O: .  ,  for  all  i,  j  =  l,...,p,  and 

For  the  study  of  the  asymptotic  non-null  distribution 
of  £,,  we  shall  consider  a  sequence  of  alternatives  converg- 
Ing  to  the  null  hypothesis  in  such  a  manner  that  the 
power  of  the  test  based  on  Z^   lies  in  the  open  interval 
(e,l).   Thus  we  replace  b^^'    in  (5-2)  by  a  sequence 
[l  +  N"^/^  p!^^3  where  p[^^  i  =  1,  .  .  .  ,p;  k  =  1,  .  .  .  ,c  are 
real  finite  constants.   Also,  we  replace  (5-5)  by 


F|^;.)(.,y)-FJJ)(.)F{5)(y) 


H(,^.)(.,y)  -H(.){.)  H(.)(y; 


+  N-l/2  4^)(.,y) 


(5.6) 


for  all  i  /^   J  =  l,...,p;  k  =   l,...,c  where  ^[j  (x,y)  are 
real  and  finite  valued  functions.   We  shall  denote  such  a 
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a  sequence  of  alternative  hypotheses  by  [h^].   Thus,  on 
defining 


GO  oo 
-00  -oo 


H(.)(x) 


J/ 


H(.)(y)  dH(.)(.)dH(.)(y) 


for  1,  j  =  l,,..,p;  k  -  l,...,c;     (5.7) 


C  =  (C:,  >  1  <  j  =  l,...,p}  ;  k  =  1, 


.c 


(5.9) 


and  following  some  routine  analysis,  we  arrive  at  the 
following  theorem. 

THEOREM  3.2.   Under  [H^  j,  £,     has  asymptotically  a  non- 
central  chl-square  distribution  with  ( c-1  )p(p+l)/2  degrees 
of  freedom  and  the  non  centrallty  parameter 


yv  -I 


k=l 


k 


'^(k)  ^*-l  ^(k)' 


(5.9) 


At    -I  JjL 

where  v  ~  Is  the  Inverse  of  v  ,  defined  In  (4.8). 

In  the  parametric  case,  the  likelihood  ratio   (L.R.) 
test  for  this  problem  Is  given  in  Anderson's  (I958)  book. 
He  has  only  considered  the  distribution  of  the  L.R.  crlterlan 
under  null  hypothesis.   However,  on  using  Wald's  (1943) 
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methods,  one  can  readily  find  out  the  asymptotic  non-null 

N 


distribution  of  the  L.  R.  crlterlan  under  {h^.   If 


Y.        ,  k  =  l,...,c  be  the  c  dispersion  matrices,  under 
{H„j,  we  will  have 

^    ~  hi       "4   ^    ~    '   k=l,...,c     (3.10) 


(k) 
where  T'  ' ,  k  =  1, . . . ,c  are  all  real  and  finite.   Now 

defining 

00   OO  00  00 


^lj,l'j'  -J     J     J     J      ^1  ^j  ^1'  ^j' 
-co  -oo  -oo  -oo 


^^l,J,l',j')(^l'^j'^l"^J')  -'^iJ^i'J'  (5.11) 

and  adopting  the  suffixing  system  of  (3.20),  we  write 

r  =  ((^rs)^r,s=l,...,p(p+l)/2     ^5.12) 
lo^^^^  (Af'    i  ^  J  =  l-'-P}.  k  =  l,...,c. 

Then,  we  have  the  following  theorem,  the  proof  of  which  is 
omitted. 

THEOREM  5.3.   Under_(5.10),  the  L.  R.  test  crlterlan 

has  asymptotically  the  non-central  chlsquare  distribution 

with  (c-l)p(p+l)/2  degrees  of  freedom,  and  the  non-centrality 

parameter 


'-f^''- 
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'  (k)  p-^  ^(k)l  (5.13) 

G    J.       O 


Comparison  of  (5-9)  and  (5-13)  will  yield  the 
asymptotic  relative  efficiency  of  the  Z^   test  with 
respect  to  the  L.  R.  test.   Evidently,  the  same  depends  on 
(i)  J.(u),  i  =  l,...,p;  (ii)  ^H^}  and  (lii)  the  unknown 
cdf  F(x) . 

Among  the  possible  members  of  the  Z-^   tests, 
the  two  most  appealing  are  (i)  the  Spearman  type  of  test 
for  which  J,^v(u)  =  /T2  (u  -  l/2),  and  (ii)  the  normal 
scores  test  for  which  J,.v(u)  is  the  expected  value  of  the 
ith  order  statistic  of  a  sample  of  size  N  from  a  standard- 
ized normal  cdf.   For  the  normal  scores  test,  it  can  be 
shown  that  if  F(x)is  also  normal,  than  (5-9)  and  (5-13) 
are  equal,  that  is,  the  two  tests  are  asymptotically 
power  equivalent.   The  detailed  study  of  the  asymptotic 
relative  efficiency  is  not  studied  here  because  the  same 
depends  on  the  unknown  F(x)  and  {Hj^I* 


6a.   NONPARAMETRIC  TESTS  FOR  hJ^^  ^ 

Under  H^^^  defined  in  (2.21)  and  (2.22),  the  basic 
permutation  argument   of  Section  5,  is  no  longer  tenable, 
and  so,  it  is  difficult  to  derive  the  permutatlonally 
distribution  free  tests.   However,  we  shall  derive  a  class 
of  rank  order  tests  by  centering  the  observations  at  the 
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respective  estimates  of  the  location  parameters  (as  Is 
generally  done  in  the  univariate  scale  problems,  see  for 
example,  Sukhatme  (1958),  Raghavachari  (I965)),  and 
working  with  these  centered  observations.   Due  to  centering, 
the  independence  is  vitiated,  and  as  a  result,  one  is  faced 
with  the  question:   Under  what  conditions  the  test  based 
on  the  centered  observations  is  asymptotic ally  equivalent 
to  the  '  'St  based  on  the  observations  centered  at  the  true 
locations?   In  this  section  we  study  this  problem.   For 
this  purpose,  we  assume  that  F^^^(x)  is  absolutely  con- 
tinuous  having  the  continuous  density  function  f^^^{x); 
k  =  l,...,c,  and  the  following  conditions  hold: 

A.  The  densities  f^^^x)  of  f(^^(x),  k  =  l,...,c  are 
diagonally  symmetric  about  their  location  parameters. 

B.  J(i)(^)  =  ^~  (")  where  f   is  symmetric  about  zero. 


C.   The  functions  4^  J , . . 

dx   (1) 


f(J!(.) 


are  bounded  for  each 


i  =  l,...,p  and  k  =  l,...,c. 

''  ( k)  {V\ 

Let  now  6^  '  be  a  consistent  estimator  of  6^  ' 

defined  in  (2.21),  and  denote 

xMk)  ^   (k)  _  g(k) 

^ia     ^ia    °i    .  1  -  i,...,P,  jg^^j 

K 

In  addition  to  the  notations  of  the  previous  sections, 
we  introduce  the  following: 
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n 


k    L- 


number   of  x\^^    such  that    {x[^^    -   b[^h 


<  X ;    a   =   1 ,  .  .  .  ,  n. 


1   =   l,...,p;    k   =   l,...,c    (6.2) 


k=l 


'N(l)^"^    -  trbr     N        ^N(i: 


(6.3) 


FfJ^^M^) 


N(1,J 


_1 
n 


k 


number   of    {x\^\    X^.J^^  )    such  that 


la 


(x(k)  _  g(k)   ^(k)  _  ^(k)) 


^  (x,y) ,  a  =  1, . .  .  ,n 


k 


.,  j=l,...,p.  {6 A) 


* 


"n(1,j)  ~  ^  ^N   ^N(i,j)   • 


(6.5) 


Now  consider  the  statistics 

OD    00 


U- 


(k) 
N(l,j) 


•^N(l 


%(l)(^l) 


'^N(j) 


!^N(j)(^j 


dF 


*(k) 
N(l,j) 


-00  -00 


1  /  j  =  l,...,p;  k  =  l,...,c   (6.6) 


m*(k) 
■^N,l 


oo 


-00 


'N(l) 


%(l)(^l) 


dFNU)(^l^^  i  =  l,...,p; 


k  =  1, .. .,c  .  (6.7) 


Finally,  concerning  5.    's,  we  shall  assume  that 
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j^l/2  i^jk)  _  g{k)|  ^^^^(1)  for  all  1  =  l,...,p; 


k  =  l,...,c.    (6.8) 
Then  the  main  theorem  of  this  section  is  the  following. 

THEOREM  6.1.   Under  the  assuinptions  A,  B  and  C,  the  random 
variables  N         LU^^ij  "  M-N^ij).  l^N,i   ^^N,i^  i  <  j  -  l,...,p 

k  =  1     cl  and  N^/2.,  (k)   _   (k)  .   ,  *(k)  _   (k), 
K  -  J.,..  .,cj  ana  iM    L^u^^^j   ^^N,ij^'  ^-^N,l    ^^N,l^' 

i  <  j  =  l,...,p,  k  =  l,...,c]  have  the  same  limiting 
distribution. 

The  proof  is  sketched  in  the  Appendix. 
It  therefore  follows  from  Theorem  6.1  that  if  we 
work  with  observations  centered  at  the  respective  estimated 

location  parameters,  the  test  considered  in  Section  4  will 

(2) 
be  asymptotically  a  distribution  free  test  for  H!^   ,  and 

will  have  the  same  asymptotic  properties. 


6b.   TESTS  FOR  hI^^ 

o 

In  an  earlier  paper  (I966),  the  authors  had  considered 

the  problem  of  testing  the  equality  of  location  parameters 

under  the  assumption  that  the  covarlance  matrices  were  the 

same.   The  tests  of  this  assumption  (that  the  covarlance 

matrices  are  same)  are  already  considered  In  the  previous 

sections.   The  purpose  of  this  section  is  to  sketch  very 
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briefly  the  problem  of  testing  simultaneously  the  Identity 
of  location  parameters  and  covarlance  matrices.   Here 
also  since  the  null  hypothesis  reduces  to  the  identity  of 
the  c  cdfs,  the  basic  permutation  arguments  of  Section  3 
will  remain  valid.   Here,  we  shall  have  to  work  with  the 
set  of  statistics  iT|.^l  .  and  t[}^'.    for  i  <  j  =  l,...,p, 
k  =  l,...,c^.   It  is  therefore  equivalent  to  work  with 
the  set  of  statistics, 

n 

y(k)       _   J^ (1)  (  j)  ,. 

^^'^^        "k  a=l   N,Rr^^  N,R\.  ^ 

for  k  =  l,...,c,  i  <  j  =  l,...,p  and 

""k  ,,,  (6.10) 


T 


k    1  ,-^  ^1    p   .   ,        1    .     ^' 
Mi  ^  TT  2_  E      ^°^  1  =  1, . . .  ,p,  k  =  1,  . .  .  ,c 

^^'^        ^\   a-1   N,R. 


The  permutation  covarlance  matrix  of  these  cp(p+3)/2 
random  variables  can  be  found  in  a  similar  manner  and  a 
test-statistic  similar  to  £,.  in  {^.2h)    can  be  constructed. 
Finally,  Theorem  5.I  may  be  used  to  study  the  asymptotic 
properties  of  this  test. 
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7a.      APPENDIX  A 


Proof   of  Theorem  5.1.      Writing 


'N( 


DrNd: 


N(J)     L^N(j)_ 


-    ^"^NCD 


H. 


J. 


N{1)J    ^N(j) 


H 


N(0). 


-   J/. 


H. 


{i)L"N(i) 


J 


(j) 


H 


N{j) 


J} 


+   J. 


(DL^Nd) 


(j 


H 


N(j). 


,     (7.1) 


(1) 


H 


N(l)_ 


J,  .. 
( j) 


H 


N(j)J 


=   J 


(1) 


H 


(1) 


J 


(j 


H 


(J)J 


■^    (^N(i)    -   ^(1))    -^(i) 


H 


(i) 


J 


(j) 


H 


(j) 


■"  ^%(j)  "  ^(j)^  ■^(j)r(j)j  '^(Dni). 


1    /  2       '  ' 


%{1)    ^    (1    -   ^^"(1), 


0H,,,  ,,    +    (1    -    6)H,, 


(j)L     N(j) 


(J) 


^|(%(j)  -«(j))' j(j)K(j)  ^(^  -^^«(J)_ 


(i) 


^Vi)    ^   ^'   -   '^"(i). 


"^^"n{i)  -  "(i)^^\(j)   -  "(J)^  J(l) 


^«N(i)    ^    (^    -'^«(i) 


0H,,,,,  +  (1  -  e)H, , 


(J)L     N(J) 


(j) 


(7.2) 
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where  o  <  0  <  1  and 


,(k) 


nj^[l,  j] 


,(k) 


dFA^'r-.   _.i  =  d(FA"r_.   _.i-  FrV'.n)  +  dF  ^  1^ 


,(k) 


n^[l,j]~   [l,j] 


(k) 

[1.  J 


(7.3) 


and,  proceeding  as  in  Bhuchongkul  (196^)  and  Purl  and 
Sen  (1966),  we  can  express 


u(k)   _   (k)        (k)        (k)       ^  (k)       ,      ^. 
"n,1J  -  ^1N(1,J)  +  ^2N(l,j)  +  ^3N(1,J)  ^   ^^rN(l,j)   ^''"^^ 

and 


,(k)  _   (k)   ,  p(k)    ,  ^(k)    ,^^(k) 


2N(1)  ^  ^^rN(l) 


(7.5) 


where 


A 


A 


(k) 
lN(l,j) 


(k) 
2N(1,J) 


00  00 

-J  J 

-00  -00 
CO  00 

^J  J 

-00  -00 


^(1) 


H(,)(x.) 


dF^^j.  .,(x.,x.)     (7.6) 
n^(i,j)^  1'  J 


%{l)^^i^  -  ^1)(^1^ 


(1) 


H(,)(x.)_ 


J 


(j 


,(k) 


hj)(^j^  ^^[i;j](^i'^j)  ^'-^^ 


>(k) 
3N{l,j) 


00   00 


-00  -00 


Kw  .,  (x  .)  -  H,  .,  (x  .) 
„N(  j)  '     j'  (  J)   J 


J,  ..   H,  .s  (x  .) 

( j)  L  ( j)   J 


J 


(1) 


H(,)(x.) 


dFf^)  .■,(x.,x.)  (7.8) 


00 


^1N(1' 


(1: 


H(,)(x.) 


f(^!,(x.)  -  pj^hx.; 

N  ( 1 )   1      ( 1 )  ^  1 


(7.9) 


-00 
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B, 


(k) 
2N(1) 


00 


L"n(i)(^i^  -  "(i)(^i) 


J 


(i: 


H(.)(x^)jur(.^V^^ 


dF!i'|(xJ(7.10) 


-co 


-A^\ 


GO  GO 


"''lN(i,j)  =]     J       K(l)  -  ^l(i)_ 


J/.N    H 


J, 


-00  -00 


(1)  r(i)i   (J) 
-(k) 


H,  .V 

.  (j1 


pv^w         (k) 
N(l,j)  -  ^(i,j).l 


(7.11) 


^: 


(k) 

'2N(l,j) 


00  oo 


L"n(j)  -  "(j 


J,  ., 
(j) 


H,  .. 


J/  . 


(j)J  "(1) 


H 


(1) 


-oo  -00 


f(^! 


,(k) 


N(i,j)  "  ^(l,j) 


(7.12) 


^(k) 
3N(1,J) 


oo  oo 


1_ 

2 


%(i)  ~  ^(i) 


2  j" 


%(1)  +  (1  -  ^)H(,)_ 


-00  -00 


(J 


,%(J) "  ^^  -  ^^^J)_ 


^%(i1j)  ^^-^^^ 


/?(k) 
'  4N(i,j) 


oo  00 


-00  -00 


H 


2  L  N(j) 


-  ^J)] 


'  j;;.)  K(  .^ .  (1  -  e)H(jj 


^(1) 


%{i)  ^  (1  -  ^)H(i)^ 


(7.1^) 


/t, 


(k) 
5N(1 


,J)  -J] 


"n(1)  -  "(i) 


"m(j)  -  "(j) 


(1) 


^N(l)^  ^N(J) 


K(i)  ^  (1  -  ^)«(i)l^;j)K(j)^(i-^)"(j) 


^^N(i,J)('^-^5) 


"  6n(i, j 
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_J(1)    ^^(1)^   J(J)    ^^J)^   ^  ^%(i)    -"(1)^ 


^^-^N(1)^^^N(J) 


^(1)  f^(i)^  ^(J)  f^(J)^  ^  (%(J)  -^(J)^ 


^(1)  f^i)^  ^(J)  t«(J)l 


dP 


(k) 
(1,J) 


(7.16) 


7y(k) 


"^Nd: 


H. 


N(1)J     ■^N(J)     L^N(j)J    '^^NCl^j)       (^-l^) 


,(k) 


^^-%(1)^%(J) 


■^Nd)  f^Nd)]  '^N(j)  f%(j)^  -  '^(i)  f%d)] 


^Nd)  n%(j 


•^(j)     ^^N(j)^ 


N(l,j) 


00 


ell^).M  =  -(N+i)-i  Tr,,,,  j). jh,, J  dpf^^ 


INd) 


'Nd)  "(i)^"(i)J  "^Nd: 


-oo 


c 


(k) 
2Nd) 


OD 


(%d)-  ^1))  j;i)f«d)^^f^Nd)-  ^d)^ 


■CO 


00 


^5Nd)    ~     j    f'^CD^NTI  ^Nd)^    '    "^(D^^d)^" 


-00 


(7.18; 


fe  ^(1)  -^d)^'^d)f^d)^^'^^N( 


1, 


(7-19) 


(7.20) 


(7.21) 


00 


'^'^Nd)  -  j  ^'^n{i)^W+i  %d)^   "  -^(D^N+T  ^(D^^^^Nd: 
-oo 


(7-22) 


Now  proceeding  precisely  as  In  Bhuchongkul  (19^4)  and 

Purl  and  Sen  (I966),  It  can  be  easily  shown  that  the 

-I/2 
'^  and  (^   terms  are  all  o  (N  '^  )  .   Thus 
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1/2  "^ 

^   '■^IN(I)  "^^2N(1)-'  t®"^s  t°  zero  In  probability 

as  N  H.  00  ,  for  each  1,J  =  1,  .  .  .,p  and  k  =  1,  .  ..,c. 
Hence  It  suffices  to  establish  the  joint  asymptotic 
normality  of 


N 


1/2 


(k) 


h= 


N 


J   ^hN(l,j)i  1  ^  1  <  J  ^  p;  k  =  l,...,c 


1/2  r^(k)    ^  ^(k)     .         ,  ,    ,      1 

L^lN(i)  +  ^2N(1)'  ^  "  l.-.-,P;  k  ^  l,...,cj 


and 


Now  let  us  denote 


la        ja 


(1) 
00  oo 


^U)<'^ 


'(i) 


."(Jj'^^Ja'']         '7-27' 


"-;L  <"' ^ // [^ixL<^.' -rM 


-co -00 


(k) 


■^(1)^^(1)^  j(j)f"(j)^  ^^(i;j)^^i'-j)    (^-^^^ 


where 


(r) 


IX^   (Xj^)  =  1  If  X.    <  X.   and  Is  zero  otherwise. 


ia 


OD  00 


u(i,j)  r^ir).    ^       r  r     r  (r)  ,   .    (r),  . 


k,r,a 


-oo  -00 


^^Ja   J  ^^J^   J 


] 


(7.2^) 


(k) 


•J(j)f"(j)^  J(i)^"(i)i  ^^(i;j) 


where 

,(r) 


(r) 


^ly   ^^1^  =  1  if  X;^  <  X.  and  is  zero  otherwise 


(7.2^) 


(7.27) 


Then,  we  can  write 
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(k) 


X 


(r) 

N 


n 


n. 


k 


^(k)   ,^(k) 
%N(1)  +^2N(1) 


N   I-  2n^,k,a^  la      3n^,k,a   ja   J 


(7-28) 


-z 


r=l 

r^k 


n 


E  B.(x!^^)}  +^  y~ 

i'  xa  -^J   n,  '^-j 
-'    k  a=l 


"(i)<4"' 


-  4"'  B.(x!l^') 

N    1   la 


(7.29) 


J 


(1) 


H,.,(x!kh1  -  A(k)B.(x!^))]} 
(1)   aa  J    N   1'  la  '  I  -> 


where 


B,(x(>^')./' 


x^(r) 


J 


(1) 


H(^)(y)   dF|jj(y) 


(7.30) 


X. 


01 


with  X  .  determined  somewhat  arbitrary,  say  by  H/.n(x  .)  =1/2; 
E  represents  the  expectation  and  X;   ,  a  =  lj,...,n,  has  the 
distribution  F.   ;  i=l,...,p;  k=l,...,c. 


i!^)  and  e[^) 
any  real  constants,  not  all  zero,  we  can  write 


Hence  with  d:.   and  e!^  ,  1  <  i  <  j  <  p,  k  =  1, . . . , c 

X  J  _L 


1/2  r^  y^  ,^^^ /„(k)       ,„(k)       ,„(k)      I 

Lltpifel      iJ      l^lN(i,j)    +^2N(i,j)    +  Vd^j)-* 
hlhl^      l^lN(i)    +  ^2N(i)iJ 


+ 


=  n1/2  ±  ,(r) 

r=l 
r/k 


n 


n 


> 


(k)  f 

u  .    .  ■) 

n     ^    I    — — :      1,1  >w     <_j.ii 

^  r  a=l    '-  i<j  k 


v-rr      ij      I.     2n,  ,r,a      la 


(x!^M  + 
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3nj^,r,a   Ja  J   -j^  1   ^-  1'  la        V    ia      \ 

^k  fcl  L  1<J  ^J  ^  l^k'«    ^ 

+  ^(k)  (k)  I   ^  (k)  r     -    ,^(k),1 


(k)>     (k)„       ,  (k). 
ja  ^  +  ^^  ^2n^,k,a^^ia  ^ 


4^)  B.(x(^^)  -  EfJ 


(1) 


(k) 


(k) 


■(k) 


"(l)'=^L'''   -C'B,(xi^^')i 


These  are  the  c-pijinmatlons .   They  Involve  Independent 
samples  of  identlCdlly  distributed  random  variables,  each 
having  finite  third  moment.   The  proof  follows. 

To  compute  the  variance  covarlance  terms,  we  note 
using  (7.24)  to  (7.29)  that 


N  var 


U 


^—   ^hN(l,j), 


h=l 


=  TZ^""^   var 
r=l 
r^k 


r2n^,r,a^^la  '    ^ 


(i.J)    ix[^))]   +   1   var  ju|i'J^x(^^x(^)) 
3nj^,r,a    ja  )        ^(k)     \  lnj^,a  ia.    '    ja 

N    2n^,k,a'  la      N   5n^,k,a   ja   / 


+ 


N  GOV  (  >   A^^^ 


A 


rj  hN(l,j)'  f^   "hN(r,s) 


-^(^iN(l,j)'  ^3S(r,s))  • 


(7.31) 


(7.52) 


(7.33) 


a=l  p=l 

The  variance -covarlance  terms  of  (n''"^^(b|}JL  ,  +  B^^L  ^ ) , 
1  =  l,...,p;  k  =  l,...,c  are  given  In  the  expression  (?.3) 
to  (5.6)  of  Purl  and  Sen  (I966).   The  other  varlance-covarlance 
terms  can  be  evaluated  analogously. 


ho 


7b.   APPENDIX  B 


The  proof  of  Theorem  6.1  rests  on  the  following 
two  lemmas  the  first  of  which  is  due  to  Raghavachari 
(1965)  and  the  second  is  a  slight  generalization  of 
Raghavachari  (I965). 

LEMMA  1.   (Raghavachari.)  If  the  assumption  (c)  holds,  then 


_d_ 
dx 


J 


H(.)(x) 


is  also  bounded,  where  H , . s ( x )  is  given  by  C^  i) 


LEMMA  2.   Let  (X^,  Y-j_ ),...,  (X^^,  Y^^)  be  i-i-d-r-vs  with 
continuous  cdf  P ( x  -  ^ ,  y  -  tj  )  where  F  is  assumed  to  be 
symmetric  about  (0,0).   Let  F^(x)  and  F^{y)    be  the 
marginal  cdfs  of  X .  and  Y .  respectively.   Let  ^^  ( X-,  ,  .  .  •  X^, ) 

A 

and  Ti  ( Y-,  , .  .  .  ,  Y^ )   be  consistent  estimators  of  ^  and  t\ 

1  /2  '^  1  /2  ^ 

respectively  such  that  N  '^  (|  -  I )  as  well  as  N  '^  (rj  -  r) ) 

is  bounded  in  probability  as  N  — >  00 .   Assume  that  there 

exists  a  function  M(x,y)  =  M^(x)M2(y)  such  that  the 

following  conditions  are  satisfied. 


(1) 


M-l(x) 


2 


(ii)  m;(x)  ^:^ 


(iii; 


<  00    E 
SM. (x) 


M2(y; 


<  00 


and  M2(y)  =  ^M^{y)/^y   exist. 


M^(x  -  t^)   <  K^T^(x),   M2(y  -   ^\<^   T2(y) 

uniformly  in  t-,    and   tp   for    |t.|<   c;    i    =   1,2;   (c .    and  K.  , 
i   =  1,2   are   constants). 


kl 


(iv)   E 


T;l(X)] 


00  ,  E 


T2(Y) 


<  00 .  where  X  and  Y 


have  the  distribution  functions  P-,(x  -  ^)  and  Fp(y  -  -n) 
respectively. 


(v)   E 
Then 


MJ(x  -  ^)j  =  0,  E  M2(y  -    r])]  = 


0. 


/   N 


Proof  of  Lemma  2.   Without  loss  of  generality  we  can 
take  ^  =  Ti  =0.   Then 


^^      N  ,       0^t.    t„   N  ,       0  t 

=  -  -i  IZ  M2(Yi)  M^(X^  -  ^  _  ^  ^  M^(^  )  M  (Y.  -  ^IJ,  ) 


(7.34) 


after  neglecting  terms  of  order  0(1/N).   Under  the 
assumptions  (1)  to  (iv),of  the  lemma,  it  is  easy  to 
note  (by  applying  Tchebychev ' s  inequality)  that 


t.  N   -       ,       e  t 
M^(Y,  )  M,  (X,  -  ^    ^ 


N 


5ru 


2'  1'   1  '^'1 


Vn 


^h 


-  e(M2(Y^)  M^  (X^  -  -i-i  )} 


v/N 


— >  0  .   (7. 3-^) 
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Furthermore,  by  dominated  convergence  theorem. 


r?   v  1 — 

e{m2(y.)  m;  (X.  -  ^J_>e|i^ 


//N 


>  EJM^(X^)  Mg  (Y^) 


(7.36) 


imlform  in  t.,    for    |t-^^|    <  c^;    and   since   under  the 
assumed  conditions,    e|m^(x^)    M^Cy^)}    =   0   it   follows   that 


sup 

hi^^i 


t^      N 


hH 


ITZ     M,(Y,)    M,    (X,    -  _ 


) 


^>0    . 


p 
Hence  the  first  factor  on  the  right  hand  side  — >Oas 

N  — >  00.   Similarly  the  second  factor  — >  0  as  N  — >  00.  Hence 
Q^  -2->  0  a  s  N  — >  00  . 

Proof  of  the  theorem.  Without  any  loss  of  generality 

we  can  assume  that  ^,  =  0,  k  =  l,...,c.   The  proof  of  the 

1/2 
theorem  is  accomplished  by  showing  that  (a)  N  ' 

(U*/!^l    -  ui^!  .)   ^>  0  as   N  — >  00    for  each  1,    j    =  l,...,p 
^    N,ij  N,ij' 

and   k  =   l,...,c,    and    (b)    ^^^^{'^l\l]    -  T^'(^l))    — >  0   a^ 
N  — >  00    for  each  1   =  l,...,p,    and  k  =  l,...,c. 

To  prove    (a),    note   that   proceeding  as   in  Theorem  1, 


we   can  rewrite  T^/ .     .  ■>    as 


*(k)    _    „*(k) 


^^(k) 


*(k) 


.f(k) 


Un,1j    =   ^N(l,j)    +  ^2N(i,j)    +  %(l,j)    +  ^'^rN{l,j) 


A 


1N(1, j) 


00     00 


-00  -00 


(1) 


H 


(1) 


(J 


H 


(j)_ 


dF 


(k) 
N(l, j) 


^3 


(k) 


00   00 


2N(l,j)   J     J 


H 


N(l)  -  "(i)   "^(1) 


H 


(1) 


(j) 


H 


dF 


(k) 


(J)r"{l,J) 


-00  -co 


3N(l,j) 


GO  00 


-OO -OO 


L"n(j)  "  "(j) 


,-1/2, 


J,  .. 
(j) 


H 


(j) 


J 


(1) 


L"(i) 


dF 


(k) 
{l,j) 


and  the  S"  terms  are  all  oJn  ^    ). 

Proceeding  again  as  in  Theorem  1,  we  obtain 


jjl/2,*(k)        *(k)      ^(k)    . 


NV2(,{k)      ^,(k)      ^  Jk) 

lN(i,j)     2N(i,j)     3N(i,j)' 


-  ^TV2 


=  N 


n 


2_  Pn   n^  ^l^U2Nkr.a^^l 
r;^k 


r  a= 


r) 
2Nkr,a^"'la 


-v(-))  + 


^5Nkr,a  ^^ja    "^j  >^         ^"2Nk,r^^la^  +  ^3Nkr,a^^ja  '^l 


+ 


n1/2  i.^  rn(i,J)(x(k)_-v  (l^J  x(^)-  v^.^^)  + 


■^k  o^T  ^  lNk,a^  la    j  ^   j 


a 


.(k)  r,j(l,J)   ,^(k)   ^(k).     (l,j)  ,^(k)   0{k)w_ 
^N   l.U2Nkk,a  ^^la  "  ^i  '    ^   "5Nkk,a^^aa   "  ^j  >} 


TT(i»J)fY(k)  y(k)x     (k)j,,(i,j)  /y(k)v  .  TT(i»  J  ) /yCk)  w  ' 
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n. 


Lp?    .  N-1/2  ^i4;'J)     (X(>)    -   v(>))    -   4J'J)     {X(^^} 

-^^    L    3Nkk,a^      ja  3         '  ?.Mlrlr    fv  ^      ir»      'J 
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B,  ^(X,  )  = 
(r) ^  ia 


4SL'^^'  '(^' 


%(x^) 


(J) 


«j(^j^ 
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00 
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J.(H.(x.))f.(x.)dx. 

J   J   J    J   J    J 


-oo 


<  K 


j|.^"[h.(x.  -  t)_ 


<  K  by  virtue  of  Lemma  1, 


Furthermore  EB,  JX   )  =  0 
I  r  ^   IQ 

Hence,  by  Lemma  2,  L^  ^  -^>  0  as  N  — >  oo 
Similarly  L^^^)^  ^{2)^^    lJJj^^  and  l|j^^  all  -^>  0  as  N 


N 


->  oo 


Hence  N 


1/2 


U 


(k) 


(k; 


N(l,j)    ^^N(l,j)' 


;  1,  j  =  l,...,p;  k  =  l,...,c 


and  N 


1/2 


u(k)        (k) 


have  the  same  limiting 


distributions.   The  proof  of  (b)  is  similar  to  that  of 
Theorem  5.1  of  Raghavachari  (I965),  and  is  therefore 
omitted. 
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